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ABSTRACT 

With its rapid response, Swift has revealed plenty of unexpected properties 
of gamma-ray bursts (GRBs). With an abundance of observations, our current 
understanding is only limited by complexity of early X-ray light curves. In this 
work, based on the public Swift data of 150 well-monitored GRBs with mea- 
sured redshifts, we find some interesting global features in the rest-frame X-ray 
light curves. The distinct spectral behaviors between the prompt emission and 
the afterglow emission implies dissimilar radiation scenarios. Interestingly, an 
unforeseen plateau is exhibited in the prompt X-ray light curves despite the 
presence of complex spikes, which might indicate the presence of a steady central 
engine. In particular, the seemingly continuous evolution with a single power law 
from the prompt to the afterglow of most GRBs might place strong constraints 
on the theoretical models. 

Subject headings: gamma rays bursts: general — radiation mechanisms: non- 
thermal 



1. Introduction 



Ever since gamma-ray bursts (GRBs) were first reported as flashes of cosmic gamma 



rays lasting for less than a tenth of a second to more than tens of seconds (IKlebesadel et al. 



19731 ). they have shown a great diversity. Unf ortunately, their energy so urce and radiation 



mechanism remain highl y speculative (e.g. , see iFishman fc Meegarull995l. for a re view) . The 



first detection of X-ray ( ICosta et al.l 119971 ) and optical ( Ivan Paradijs et al.l 119971 ) afterglows 
that faded as a sing le power law oc strongly favored the prior relativistic fireball models 
( jWijers et al.lll997l ). Even though the origin of prompt gamma-rays remained uncertain, the 
longer-wavelength afterglows, especiall y the optical af terglows, have been found broadly 



consistent with fireball models (e.g., see iMeszarosI |2002| . for a review 
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The fundamental under standing; of GRB afterglows has been questioned by the Swift 
mission (iGehrels et al.1 12004). w hich has revealed several unexpected behaviors in the X-ray 



afterglows (e.g., see IZhang||2007l . for a review). Some bursts have an initial steep de cline, fol- 



lowed by a shallow decline a few hundred seconds late r (e.g., Tagliaferri et al. 20051). Others 



have erratic flares with strong spectral variations (e.g.. Burrows et al. 2005];^ 



20071). As a group, they were proposed to have some canonical behaviors (jZhang et al. 



Nousek et al. 



Chincarini et al. 



2006 



2006 



power- 



flCostalll999l ). 



Contradictorily , some other bursts also show evidences for a single 



law d ecline (ILiang et al.l 120091 ). which are consistent with the pie- Swift observations 



With the X-ray afterglows being puzzlingly diverse, it is a great challeri ge to produce 



an applicable and self-consistent physical understanding (e.g., see lZhang||2007l . for a review 



P revious works on GR B diversity were based on relatively restricte d samples, e.g., 27 bursts 



m iNousek et al.l ( l2006l ). with only 10 measured re dshifts; 33 bursts in IChincarini et al.l ( 120071 ). 
with only 9 measured redshifts; and 19 bursts in lLiang et al.l (120091 ). with only 12 measured 
redshifts. With more than 500 Swift GRBs detected up to now, we realize that it is urgent 
to revisit the diversity issue with an extended sample. 

Using a collection of 150 Swift GRBs with well detected X-ray afterglows and known 
redshifts, we find that there exist some underlying global features in the rest-frame X-ray 
light curves, which might clarify the diversity issue and provide some strong constraints 
on the theoretical models. The structure of this Letter is as follows: the sample and data 
analysis are presented in Section [21 the rest frame light curves are interpreted in Section [31 
and the discussions and conclusion are given in Section [H 



2. Sample and Methodology 



Swift is a multiwavelength observatory (IGehrels et al.l 120041 ). which carries three in- 
struments: Burst Alert Telescope (BAT), X-Ray Telescope (XRT), and Ultraviolet /Optical 
Telescope (UVOT). We made extensiv e use of the autom ated BAT- XRT products provided 
by the UK Swift Science Data Centre (lEvans et al.ll2010l ). the online BAT GRB Event Data 
Processing Report provided by Taka Sakamoto and Scott D. Barthelmy, and the online big 
table of all well-localized GRBs provided by Jochen Greinei0. 

In Table [H we hst the 150 bursts from GRB 050126 to GRB 100425A with measured 
redshifts. We find the mean redshift < z >~ 2.14 (with the median ~ 1.93), which is 



http://www.inpe.nipg.de/^jcg/grbgen.html, which is mainly our source for measured redshifts. 
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significantly smaller t han the previously evaluated (e.g., < z >~ 2.8 for 16 bursts by 
Jakobsson et al.ll2006[ ). The distributions of the BAT duration Tgg both in the observer 
and rest frame are shown in Figure [TJ The mean value o f Tgn in observer ' s frara e is 80.4 s, 
consistent with that deduced in a sample of 237 bursts by lSakamoto et al.l ( 120081 ). With the 
150 measured redshifts in our sample, we find that the mean value of Tgo in the rest frame is 
29.2 s. Shorter bursts increase significantly in number, but no positive classification can be 
made. Even so, as indicated in Figure [H there is an apparent oddball, GRB050509B, which 
has an extremely short duration Tgo ~ 0.04 s in the rest frame. 

To view the overall behaviors of GRBs in their rest frame, we take advantage of the 
Swift Burst Analyser, which provides the combined B AT-XRT light cur ves in the form of fiux 
density Fy{v) at 10 keV in the observer's rest frame (lEvans et al.ll2010[ ). With the measured 
redshift the isotropic spectral luminosity at a given photon frequency vq (which is 10 keV 
in this work) in the rest frame can be calculated by 



i-r(t) 



AtiF,{uoM)DI{z){1 + 



,r(to)-2 



1^ 



where z/g = + z) is the emitted photon energy in the rest frame for the observed 
photon energy and t = to/{l + z) is the time measured in the rest frame related to 
the time to measured in the observer's frame. Here, a single power-law spectrum has been 
assumed F^(u, t) oc u^~^^^\ whe re T(t) is the time- dependent photon index, which is available 
by the Swift Burst Analyser ( lEvans et al.l 12010) . Throughout this work, the luminosity 
distance -Dl(^) is calculated by assuming the cosmological parameters Hq = 71 km s~^ Mpc^^, 
Qm = 0.27, and Q\ = 0.73. According to the error propagation rule, the uncertainty in the 
isotropic spectral luminosity can be given by 



(2) 



where ap^ and ur are the uncertainties in the measured fiux density and photon index, 
respectively, and the uncertainty in GRB redshift measurements is considered negligible. 



3. Underlying Global Features 

In FigureEl the calculated X-ray light curves from 0.01 s to 10'' s after BAT trigger at 10 
keV in the rest frame with the spectral evolution of 150 Swift GRBs are plotted together. In- 
terestingly, some underlying global features are revealed in both the light curves and spectral 
evolution, despite the fact that some bursts may have some non-trivial behaviors. The bulk of 
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the light curves build up a plateau lasting for about tens of seconds with significantly spectral 
softening and a following power-law decline with almost invariant spectra until the detection 
limit is reached. This feature, for instanc e, is well exhibited by the remarkable naked-eye 
burst GRB080319B@ ( Racusin et al. 2008 ). There are also several apparent outliers, e.g., 
GRB060218, GRB100316 D, and GRB060614. The former two showed clear evidences of an 



associated supe rnova (e.g.. ISoderberg et al.ll2006l : IChornock et al.ll2010l : IStarling et al.ll2010 



Fan et al.l 120101 ) , while o n the contrary the latter showed unexpec ted evidences of no as- 
sociated supernova (e.g., iFynbo et al.ll2006l : iDella Valle et al.ll2006l ). Another two peculiar 
ones are the extremely short burst GRB050509B and the optically dark burst GRB090417B. 
The former had an overall much lower luminosity, and the latter s howed evidences of d ust 
extinction a nd scattering , revea led b y the strong spect ral softening (IHoUand et al.ll2010l ). as 
predicted by lShao &: Pail (120071 ) and lShao et al.l (120081 ). Interestingly, those outliers tend to 
have a lower redshift. 

To explore the underlying global features mentioned above, we fit a broken power- 
law model to the light curves with the 145 GRBs (i.e., without the apparent outliers, 
GRB050509B, GRB060218, GRB060614, GRB090417B, and GRB100316D), assuming a 
four-parameter function 



t 

*L,b 



«1 



02 



(3) 



where {tL,h,Lu^b) is the transition position between the two power laws with the indices as ai 
and a2, respectively. The best-fitting results (top panel in Figure [31 dashed hne) are t^^b = 
15.8±0.4 s, L^,b = (1.2±0.1) X 10=^^ g-i^^-i, = -0.04±0.01, and aa = -1.47±0.01. 
The goodness of fit is measured by the adjusted coefficient of determination, ~ 0.99928, 
which is very close to 1.0 and indicates a good fit. Separately, we fit the spectral evolution 
with another four-parameter function 



Ffe + 5i log (^^) , t<tr,b, 
Fb + 52 log f T^- V t>tr,b, 



(4) 



where (tr.bjTb) is the transition position between two linear segments in the log-linear plot, 
with the slopes as 6i and 62, respectively. The best fitting results (bottom panel in Figure |3l 



^The previously reported discontinuities between the BAT-XRT combined light curves have been min- 
imized by a c count ing for the spectral evolution. For details of how these light curves were produced, see 
Evans et al.) |2010l ). 
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dashed line) are tr,b = 131.2±4.2 s, Tb = 2.10±0.01, Si = 0.47±0.01, and 62 = -0.03±0.01 
~ 0.93809). 



At this point, based on the hght curves and spectral evolution history, two emission 
ingredients can be distinguished in the X-ray emission: an early plateau with significant 
spectral softening and a following power-law decline without spectral variance. Given that 
tL,h ~< ^90 >; this strongly suggests that the former is connected with the prompt emission 
and the latter is connected with the afterglow emission. Since these two components have 
been thought to arise from different radiation regions in the relativistic ejecta, the seemingly 
nice transition in the light curves may raise a fine-tuning issue. Nevertheless, the discrepancy 
between t^ b and tr,b implies that the interpretation is more complex. 

To further explore the details at the junction, we now introduce two additional breaks 
both in the light curve and in the spectral evolution. Therefore, we have an eight-parameter 
model for the light curve defined by 



where, (ti,Li) and (t^jL^) are the positions of the first and third break points, respectively, 
and ai, 0^2, 0:3, and are the sequential four power-law indices. Given these eight free 
parameters, the position of the second break point is fixed according to Equation ([5]) and 
can be determined by 



The best-fitting results (top panel in Figure [3l thick line) are ti = 26.5 ± 0.4 s, Li = 
(1.4 ± 0.1) X 10^2 erg s-^Hz-\ = 620.3 ± 80.9 s, L3 = (4.3 ± 0.7) x 1029ergs-^Hz-\ 
ai = 0.02±0.01, a2 = -2.58±0.03, as = -0.90±0.07, and 04 = -1.35±0.01 {R^ ~ 0.99932). 
Given the best-fitting parameters, we have t2 ~ 154.6 s, and L2 ~ 1.5 x 10^° erg s~^Hz~^ 
according to Equations and ([7]). 





(6) 




(7) 
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Similarly, we have another eight-parameter model for the spectral evolution defined by 



where (ti,ri) and (tsjFs) are the first and third break points, respectively, between linear 
segments in the log-linear plot, and 6i, 62, 63, and 64 are the sequential four slopes. Given 
these eight free parameters, the position of the second break point is also fixed according to 
Equation ([8]) and can be determined by 



The best-fitting results (bottom panel in Figure |3l thick line) are ti = 13.4 ± 0.7 s, Fi = 
1.45 ± 0.01, ts = 421.8 ± 29.2 s, F3 = 1.99 ± 0.01, 61 = 0.17 ± 0.01, 62 = 0.79 ± 0.01, 
^3 = -0.47 ± 0.05, and ^4 = 0.04 ± 0.01 (i?^ ~ 0.93993). Given the best-fitting parameters, 
we have ^2 ~ 131.0 s and F2 ~ 2.23 according to Equations and (ITU]) . 

The goodness of fit has been improved for both the light curve and the spectral evolution, 
though the spectral evolution still appears to be a little more deviant than is presumed by 
the fitting model. The general agreement of the three break times of the light curve and 
those of the spectral evolution, in spite of the fact that the two fittings are made separately, 
justifies the existence of the extra breaks we have just assumed. Therefore, the four-segment 
fitting model (with eight free parameters) is a more realistic interpretation for the underlying 
global feature, which, hereby, is composed of an early plateau (oc t°) with a mild softening 
(AF ~ 0.5 on average), a following steep decline (oc t~^-^) with a further severe softening 
(AF ~ 1.0 on average), a newly emerging shallow decline (oc t~^-^) with a slight spectral 
hardening (AF < 0.5 on average), and a following single power-law decline (oc t~^'^) without 
spectral variation (F ~ 2) until the detection limit is reached. However, as indicated by 
Figure [3l the four-segment model is not overwhelmingly better than the two-segment model, 
since the steep decline and the shallow decline of most bursts are interlacing with each other 
and cannot be well defined in both the light curves and spectral evolution. For some bursts 
(e.g., GRB080319B), the existence of the shallow decline and the steep decline is debatable. 
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(10) 



(9) 
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Discussion and Conclusion 



GRBs have been exposed to be diverse botli in their observed prompt emission and 
afterglow emission. As a result, the physical interpreta tion has remained controversial. In 
this work, based on the X-ray data (lEvans et al.ll2010l ) of 150 well-monitored Swift GRBs 
with measured redshifts, we find, in the rest frame, some underlying global features concealed 
in the observed complex X-ray light curves and the spectral evolution, which could clarify 
the diversity issue. 

Generally speaking, two distinct and consequent emission ingredients unequivocally exist 
in the X-ray light curves as explicitly revealed by the spectral differences shown in Figure 
[3l which can be simply recognized as the prompt emission with apparent spectral softening 
and the afterglow emission without apparent spectral evolution, respectively. The most 
straightforward global features are the early plateau (oc t"^ with ai ~ on average) with 
a mild spectral softening (with a photon index increase AF ~ 0.5 on average) and the late 
single power-law decline (oc t"" with 0:4 ~ —1.4 on average) without spectral variation (with 
a constant photon index F ~ 2 on average) up to t ~ 10'' s, which are the previously known 
prompt emission and afterglow emission, respectively. 

The transition from the prompt emission to the afterglow emission undergoes two dis- 
tinct phases: a steep decline (oc t"^ with 02 ~ —2.6 on average) with a severe spectral 
softening (AF ~ 1.0 on average) and a shallow decline (oc t°"^ with 0:3 ~ —0.9 on average) 
with a mild spectral hardening to agree with the invariant spectra in the afterglow emis- 
sion. The temporal feature is consiste nt with the canonical behavior previously proposed 
(IZhang et al.ll2006l : iNousek et al.ll2006l ). The spectral feature supports the speculation that 
the steep decline is connected with the prompt emission, while the shallow decline is more 
relevant to the afterglow emission. On the other hand, the erratic X-ray flares reportedly 
existing in half of the bursts are not self-explanatory in our analysis. Since they exhibit 
significant spectral softe ning and their peak lu minosity follows a steep power-law decline 
Lpk oc statistically ( IChincarini et al.ll2010l ). they should be connected with the steep 
declines, which also explains their analogy to the prompt emission. 

Furthermore, as proposed earlier the spectral feature could be a crucial signature for 
discriminating the prompt emission from the afterglow emission. For this reason, the high- 
energy emission detected by the Fermi satellite which normally decays as a single power 
law with an index a ~ —1.4 and shows no significant spectral evolution could be more 



convincingly interpreted as the emission of the external forward shock (e.g.. lZou et al.ll2009 



Kumar fc Barniol Duran II2009I : iGao et al.ll2009l : iGhisellini et all 120101 ). 



Interestingly, despite being composed of complex spikes, the superposition of the prompt 
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X-ray light curves of different bursts exhibits an unexpected plateau lasting about ti ~ 30 s 
on average. This might signal a rather steady central engines. Furthermore, it indicates that 
the variabilities displayed by the complex spikes may just manifest the internal instabilities 
in the relativistic outflows launched by the central engines. For the family of black-hole 
central engines, this may shed some lights on the process of accreting onto a nascent black 
hole. We assume that the accretion process does not significantly alter the mass of the 
black hole (Mbh) and the accretion timescale is that of the free fall t ~ i?^/^A/GMbh ~ 
50 s -R^o^(Mbh/2M0)~^/^, where G is the Newton's constant and R is the initial distance of 
the accreting matter to the central engine. Then we simply estimate the accretion rate as 
dM/dt oc nR^dR/dt oc rai?^/^ where n is the number density of the star material being 
accreted. The "steady" prompt emission may suggest that dM/dt oc t*^, which could in 
turn imply that n oc R~^^'^. Such a density profile should hold up to a radius i? ~ 5 x 
10^ cm (Mbh/2M0)^/^, beyond which the density profile should be steepened to n oc R~^ '^ 
to account for the steep decline oc supposing the luminosity of the X-ray emission is 

still proportional dM/dt. The steep decline continues up to t2 ~ 155 s corresponding to 
i? ~ 2 X 10^*^ cm. On the other hand, for a magnetar-like central engine, the plateau phase 
can be naturally accounted for if the prompt emission is powered by the dipole radiation of 
the r aagnetar L^m (x (l-|-t/rn)~^ oc t° for t -C tq, where tq is the so-called spin-down timescale 



(e.g.. lPacinilll967l : lGunn &: Ostrikerlll969[ ). The quick decline phase could indicate that either 
the magnetar has collapsed due to significant accretion or the efficiency of converting the 
wind energy of the spin-down {t > tq) magnetar decreases with time. 

Last but not least, we need to mention that the seemingly continuous evolution with 
a single power law from the promp t to the afte rglow phase of most GRBs, which was first 
hinted by the BeppoSAX satellite (jCostal 119991 ). might provide some strong constraints on 
the theoretical models. Otherwise, a severe fine-tuning would be needed. Meanwhile, the 
well-known canonical behavior seems to be overemphasized. The presence of a steep decline 
followed by a shallow decline (oc with a > —0.5) in the early X-ray light curves might not 
be as prevailing as previously thought. Even though, a statistical bias may exist with respect 
to the bursts with a well-determined re dshift, which should be brighter and less affected by 
dust extinction (e.g.. lFynbo et al.ll2009l ). Interestingly, the bursts in our sample that exhibit 



a single power-law decline tend to have a relatively larger luminosity. While the bursts that 
have puzzlingly different features, especially those that clearly exhibit a steep decline followed 



•^In this scenario, the steep X-ray dechn e is attributed to the weakening of the activity around the central 
engine, as propos ed bv iFan fc Weil (|2005[ ). In the high-latitude emission (i.e., curvature effect) model (see 



Zhang et al.ll2007i and references therein), the decline should be (t — to) , where to is the ejection time 



of the last main pulse in which high latitude emission dominates over that of all other pulses. For the global 
feature of the steep decline found in Figure |3] this model works for to ~ 0. 
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by a shallow decay, tend to have a much lower luminosity during the shallow decline, which 
may give some hint on the origin of the shallow declines and disfavor the energy injection 
models. 

We are very grateful to the referee for insightful comments and Z.-G. Dai, B. Zhang, E.- 
W. Liang, S. Covino, and Z.-P. Jin for helpful communications. L.S. is grateful to N. Mirabal 
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Table 1. List of 150 Swift GRBs with Measured Redshift. 
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Fig. 1. — Distribution of the BAT burst durations Tgo measured in observer's (dashed his- 
togram) and rest (sohd histogram) frame, respectively. 
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Fig. 2. — X-ray light curves (top panel) and spectral indices (bottom panel) of 150 Swift 
GRBs in the rest frame. Six individual bursts with different redshifts are indicate d with 
different colors from top to bottom. The rest are plotted in gray color (jCapridl2005l ). 
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Fig. 3. — X-ray light curves (top panel) and spectral indices (bottom panel) of 145 Swift 
GRBs (described in the text) in the rest frame. The best-fitting models with four parameters 
and eight parameters are shown with dashed lines and solid lines, respectively. 



